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Abstract This study was carried out to investigate the 
postbuckling behavior of thin wood-based sandwich panels 
under high humidity. Using the Rayleigh-Ritz method 
based on the von Karman nonlinear theory for the panel, 
the solutions for both the approximate and the closed 
form for postbuckling of orthotropic panels were derived to 
evaluate the deflection for the boundary condition of all 
clamped edges. The results suggested that the edge move¬ 
ment be considered for evaluation of a critical moisture 
content and deflection of thin wood-based panels fixed on 
the core with an adhesive. The numerical solution obtained 
from the derived model showed some discrepancy with the 
experimental results. The predicted results overestimated 
the center deflection of the panels because creep and plastic 
deformation might be caused by considerable in-plane 
stress on panels. 
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Introduction 

Most thin wood-based panels used for door skins, walls, and 
ceilings are attached to the core materials with nails or 
adhesives and are used for interior building applications 
and outdoor applications where the temperature and hu¬ 
midity fluctuate. The dimensional changes of these wood- 
based panels are likely due to moisture content (MC) 
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changes, which can go up a few percentage points within a 
day . 1 At a critical MC, the panel loses its dimensional stabil¬ 
ity, which leads to hygrobuckling. The surface waviness of 
wood-based panels may be due to out-of-plane displace¬ 
ment after hygrobuckling begins. This problem occurs often 
with the decorative panels coated with high glossy paint, 
which damages its decoration. Thus, it is necessary to un¬ 
derstand this phenomenon on the basis of postbuckling 
theory. 

Numerous studies on the buckling of laminated com¬ 
posite plates have been done, and Leissa 2 has provided an 
excellent review on the subject. For example, the 
postbuckling of anisotropic plates subjected to mechanical 
forces was studied by the Airy stress function and beam 
characteristic function in conjunction with both clamped 
and simply-supported boundary conditions . 3 A closed-form 
solution for simply-supported conditions and a finite ele¬ 
ment method for simply-supported and clamped conditions 
were reported for beam conditions . 4 Compared to those 
on mechanical forces, however, there are few studies on 
the thermal postbuckling behavior of laminated plates . 5,6 
The closed-form approximate solution for orthotropic 
square plates with simply-supported edges were reported 
using the Rayleigh-Ritz method . 7 Studies on thermal 
postbuckling have made gradual progress in recent decades, 
but most of the studies focused on plates with simply- 
supported edges .^ 12 

From the literature review, the earlier studies on thermal 
postbuckling behavior could provide a tool for better un¬ 
derstanding the postbuckling behavior of orthotropic 
plates. However, the earlier studies could not be used to 
predict postbuckling of orthotropic plates because the 
boundary conditions of wood-based sandwich panels are 
assumed to be clamped edges, and these panels also have 
movable edges owing to the expansion of facing panels. 

To the authors’ knowledge, the closed-form solution for 
orthotropic plates with clamped edges is not available. At 
the same time, few attempts have been made to find ap¬ 
proximate solutions for columns. Moreover, some research 
has been done on the postbuckling of wood-based panels, 
such as that on strips . 11 
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To optimize the structure of panel application, including 
the distance between the core material and surrounding 
conditions, this paper examined the hygropostbuckling be¬ 
havior due to hygroexpansion of hollow-core wood-based 
sandwich panels composed of hardboards for the faces and 
laminated veneer lumber (LVL) for the core. The relation 
of hygroexpansion and stress in wood-based material is not 
elastic but is viscoelastic and mechanosorptive under mois¬ 
ture movement. This study was based on elastic behavior 
for the basic research. 

For this purpose, both double sine series of the dis¬ 
placement function and the Rayleigh-Ritz method were 
employed to develop numerical models. The derived ap¬ 
proximate solutions were compared to the experimental 
results of sandwich panels. The closed-form approximate 
solution for orthotropic panels was derived using only one 
term of a series for a reference. 


Theory 

The total potential energy of an ortho tropic panel 13 sub¬ 
jected to biaxial hygroexpansion forces is given as 
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where 

N? = AMC(a x A n + a y A n ) 
N” = AMC(a x A 21 + a y A 22 ) 
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See Appendix for explanation of the abbreviations and 
symbols. 

Nonlinear strains in midplane of panels defined by von 
Karman theory are given as 
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The stress-strain-moisture content relations for 
orthotropic panels is expressed as 
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There are four boundary conditions for all edges of the 
clamped panel, 14 but one was selected as shown as Eq. (5) 
for this work because it was consistent with immovable 
edges that were rigidly clamped to the panel. 

w = 0, u = 0, v = 0, dw/dx = 0 for x = 0, x = a ^ 5 ) 

w = (),« = 0, v = 0, dw/dy = 0 for y = 0, y = a 

These assumed boundary conditions should be noted for 
comparison with other results. Deflection equations satisfy¬ 
ing the above boundary condition in each direction are 
expressed as 15 
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The displacement functions can be obtained for an infi¬ 
nite strip, v = 0, by using Eq. (7) derived from Eq. (6) 
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Substituting Eqs. (6) and (7) into Eq. (1), performing the 
integration, and minimizing the energy expression leads to 
Eq.(8). 

-^- = 0, = o, |iL = o (8) 

dA mn dB mn dC mn 

Equation (8) gives simultaneous algebraic equations, 
which can be solved to obtain the Ritz coefficients, A m , B mn , 
r 

mn• 


Experiments 

For construction of a test assembly, a hardboard panel with 
a density of 0.96 g/cm 3 was selected as facing material be¬ 
cause it has been mainly used for door skins in industrial 
practice. The length/thickness ratios ( a!h ) were adjusted to 
200:3, 300:3, and 400:3, respectively. The dimensions of 
the prepared panel (a X b) were 200 X 200 mm, 300 X 
300 mm, and 400 X 400 mm with a constant thickness of 
about 3 mm. Specimens for determining the modulus of 
elasticity (MOE), dimensions, and moisture content (MC) 
changes were cut along the four edges of the buckling speci¬ 
men. Linear expansion (LE), for which the span was 
200 mm, was measured with digital Vernier calipers 
(±0.01 mm). The core material was laminated veneer lum¬ 
ber (LVL) made of radiata pine with a dimension of 26 X 
30mm (width X thickness). All the specimens cut from the 
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Fig. 1. Geometry of sandwich panel and assembly setup for test (units 
are millimeters) 


panel were preconditioned at 25°C and 60% relative 
humidity (RH), which is below the indoor equilibrium 
moisture content (EMC). However, the core LVL was kept 
indoors and coated with oil-based paint to minimize the 
effect of the dimensional change on the test. As shown in 
Fig. 1, the square sandwich panels were symmetrically 
assembled with polyvinyl acetate (PVA) adhesive and 
were conditioned in a constant-humidity chamber of 25°C 
and 90% RH. Holes with an 8 mm diameter was made at 
the middle of the core lumber to reduce the difference in 
the equilibrium rate between sample panels for buckling 
and dimensional change specimens. Strip specimens were 
fixed on the steel beam with a nut and bolt, which was 
enough to resist hygroexpansion forces. The center 
deflections of the sandwich panel and strip subjected to 
increasing relative humidity were measured at some inter¬ 
vals by a dial gauge (Fig. 1) and were compared with the 
numerical and analytical solutions derived. For the analysis, 
it was assumed that the MOE and the thickness change 
linearly with MC. 


Results and discussion 

Changes of MC and material properties 

Figure 2 shows the changes of MC and LE as a function of 
time at 90% RH. It took about 130 h for hardboard to reach 
an equilibrium state. Depending on the type of span used, 
the control specimens showed some variations in the 
changes of MC, thickness, linear expansion (LE), and MOE 
because they were cut from different panels (Table 1). The 
average changes of MC, thickness, and MOE were 4.6%, 
4.3%, and 11.9%, respectively. The coefficients of LE for 
different span types ranged from 0.397 X 10 -3 to 0.461 X 
10~ 3 (average 0.425 X 10' 3 ). 



Fig. 2. Change of moisture content (MC) and linear expansion (LE) of 
hardboard with time at relation humidities (RHs) of 60%-90%. Solid 
lines , /1MC; dotted lines , LE; circles, a = 200 mm; squares, a = 300 mm, 
triangles, a = 400 mm 



0 12 3 4 

W c /h 

Fig. 3. Nondimensional load-center deflection curves of clamped 
square panels with immovable edges 


Nondimensional center deflection of square panels with 
immovable edges 

The first four terms in each of the series for u, v, and w were 
used for the numerical analysis. Equation (8) gave a set of 
12 simultaneous nonlinear equations for A mn , B mn , and C mn . 
To solve these equations, Newton’s method was adopted 
because it always converges if the initial approximation is 
sufficiently close to the solution. It converged within 10 
iterations for the present analyses, giving a convergence 
criterion of 1 X 10“ 5 . The values for the nondimensional 
center deflection, W c , of square plywood and hardboard 
were obtained using the first four terms and were compared 
with those obtained using the first term, as shown in Fig. 3. 
It was assumed that hardboard was an isotropic, homoge¬ 
neous panel with a Poisson’s ratio of 0.25. For the plywood, 
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Table 1. Physical and mechanical properties of hardboard at 25°C, 60% RH 


Mechanical properties 

Span 200 mm 

Span 300 mm 

Span 400 mm 

Moisture content (%) 

6.79 ± 0.58 

6.73 ± 0.30 

6.46 ± 0.48 

(10.68 ± 0.71) 

(11.25 ± 0.69) 

(11.81 ± 0.84) 

Density (g/cm 3 ) 

0.97 ± 0.02 

0.96 ± 0.02 

0.95 ± 0.01 

Thickness (mm) 

2.94 ± 0.01 

2.91 ± 0.01 

2.92 ± 0.03 

(3.05 ± 0.02) 

(3.02 ± 0.04) 

(3.08 ± 0.04) 

MOE (GPa) 

4.57 ± 0.53 

4.22 ± 0.74 

4.38 ± 0.58 


(4.02 ± 0.33) 

(3.72 ± 0.68) 

(3.86 ± 0.51) 

Coefficient of LE (Xl(T 3 ) 

0.461 ± 0.047 

0.397 ± 0.096 

0.418 ± 0.065 


Results are averages of 21 to 28 tests 

Numbers in parentheses represent the value at 25°C and 90% RH 
RH, relative humidity; MUE, modulus of elasticity; LE, linear expansion 


which is one of the typical orthotropic wood-based panels, 
the assumed E X IE 2 , G 12 IE 2 , and v u values were 5.0, 0.2, and 
0.12, respectively. 

For the square panel, the coefficient C u was dominant, 
and A n was dominant among A mn . For hardboard, A m „ was 
equal to B nm . Figure 3 showed that the center deflection 
decreased with increasing MOE. As shown in Fig. 4, maxi¬ 
mum tension strain occurred at the midpoint (all and bll) 
and the compression strain at all and b/8 or a/2 and lb 18. In 
addition, the one-term displacement function was chosen to 
develop the closed-form solution, shown as Eq. (9). 
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Adopting the same procedures, the solution resulted in 
Eq. (10) and is shown in Fig. 3. 
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Table 2. Comparision of nondimensional center deflection with that of 
Chia for square plates 


NJN a 

Hardboard 


Plywood 


Present 

(4.661) 

Chia 

(3.831) 

Present 

(10.024) 

Chia 

(9.407) 

1.1 

0.438 

0.659 

0.421 

0.789 

1.5 

0.976 

1.546 

0.929 

2.117 

2.0 

1.376 

2.329 

1.293 

3.529 

2.5 

1.681 

3.152 

1.556 

4.879 

3.0 

1.939 

4.161 

1.759 

6.327 


Numbers in parentheses represent the value of the nondimensional 
critical load, N„ 


y 



Despite of the differences in the displacement functions, 
the values of the center deflections obtained using Eq. (10) 
were almost the same as the one-term solution from Eq. (6). 
Also, the N a values obtained from both Eq. (10) and 
the four-term Eq. (6) were 4.679 and 4.661 for isotropic 
square panels, respectively. However, the differences were 
small. 

To determine the influence of using different boundary 
conditions for the prediction of N a and W c , the present 
solutions were compared with Chia’s nine-term solutions 3 
which were obtained with boundary conditions different 
from the one used for this work. Chia assumed that the 
loaded and unloaded edges were uniformly displaced by 
normal stress. Table 2 shows that the present solutions for 


Fig. 4. Typical nonlinear strain, ty, of panels 


N CI are slightly larger than Chia’s results, but the VE C values 
were much less than those of Chia. 

Postbuckling of the hollow core sandwich panel 

Figure 5 shows the experimental and predicted results of 
the postbuckling of the hollow core type sandwich panel. If 
was not clear whether test specimens for the sandwich 
panels were flat during preparation because it was difficult 
to measure the extent of the initial flatness throughout the 
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Fig. 5. Comparison of the predicted and observed center deflections of a sandwich panel with movable edges 


panel. Therefore, some errors due to the presence of imper¬ 
fections might have influenced the predicted results. This 
was evident in the case of the short span (a = 200). 

The evaluation of W c using the nondimensional center 
deflection for panels with immovable edges is shown in Fig. 
3 and was much larger than the observed results. This might 
be attributed to the in-plane stress whose distribution was 
uniform in-plane before initial buckling and nonuniform 
beyond buckling, as shown in Fig. 4. Therefore, its edges 
might have movements due to hygroexpansion forces on the 
surface of panel, and the edge movement should be taken 
into consideration for the prediction to accommodate 
errors. 

The edge movement due to the LE of facing panels can 
be calculated by modifying the theory of solid laminates. 16 If 
the LE of the core is ignored and edges are displaced uni¬ 
formly, it can be represented as Eq. (11). 
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where the average strains in each direction are expressed as 
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Therefore, the effective hygroexpansion forces are given 


as 


N x = (a x AMC — e, )A X ( + [a y AMC — e 2 )A I2 

Ny = (a x AMC - £, )yf 2I + [a y AMC - e 2 )A 22 (13) 

The average strain that contributed to the edge move¬ 
ment was represented as a function of LE in Fig. 6 in which 
the MOEs parallel to the grain and perpendicular to the 
grain of the LVL core were assumed to be 10.0 and 0.5 Gpa, 
respectively. As expected, the average strain was inversely 
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Fig. 6. Predicted maximum and minimum membrane stress (o x ), and 
ratio of average strain to LE of sandwich panels with movable edges. 
Solid lines, membrane stress; dotted lines, ratio of average strain to LE; 
circles, a = 200 mm; squares, a = 300 mm, triangles, a = 400 mm 


proportional to the span and decreased as deflections in¬ 
creased. The center deflections obtained by considering the 
edge movement by Eq. (13) were closer to the average 
observed than those by Eq. (1), but all the predicted results 
overestimated them. The reason might be the plastic defor¬ 
mation at high MC and high stress level, 18 as shown in Fig. 6. 
It was clear from the results that the errors increased as the 
span increased because membrane stress was inversely pro¬ 
portional to the edge movement. In the cases of a = 300 and 
400 mm, it was estimated that tensile stress occurred at both 
points a/2 and b/2. 


Postbuckling of a strip with immovable edges 

The first five terms in each of the series of Eq. (7) were used 
to predict the center deflection and membrane stress. Simi¬ 
larly, Newton’s method was used to solve a set of 10 simul¬ 
taneous nonlinear equations for the panels. The analysis 
showed that the coefficient C 1 was most dominant, and A 2 
was dominant among A m . The rest of the coefficients, except 
C x and A 2 . were close to zero. 

Buckling stresses for span lengths of 200, 300, and 
380mm were 3.5,1.4, and 0.91 MPa, respectively. The stress 
increased as the hygroexpansion increased in the panel, but 
there was little change beyond the buckling of the strips. 

The center deflections predicted by Eqs. (7) and (8) were 
independent of the MOE, which was different from that of 
sandwich panels. It was almost the same as those of the 
measured deflection, as shown in Fig. 6, particularly for 
span lengths 300 and 400 mm. However, with a span length 
of 200 mm the predicted deflection was larger than the mea¬ 
sured result, which could be attributed to the plastic defor¬ 
mation of a large buckling stress and the center deflection 
beyond buckling. 

The center deflection due to hygroexpansion can be also 
calculated. 417 To understand the postbuckling behavior, an¬ 
other equation can be derived using the relation between R 
and e, as shown in Eq. (14); 



Fig. 7. Comparison of the predicted and observed center deflections of 
a strip. Solid lines, Eq. (16); dotted lines, Eqs. (7) and (8): solid lines 
and symbols are averages; circles, a = 200 mm; squares, a = 300 mm, 
triangles, a = 400 mm 
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Assuming that once a critical load is reached and the strip is 
incapable of supporting any further load, the strain is given 
as: 


rr 2 h 2 

s = a x AMC-——y- (15) 

3 a~ 

Therefore, the center deflection can be calculated by replac¬ 
ing the R value given by Eqs. (14) and (15) with Eq. (16). 
The calculated value is about 4% smaller than Spalt’s 
result. 1 ' 


v x=a! 2 


8 R 


(16) 


The center deflections obtained by the above equations 
are represented in Fig. 7. The results show the same trend as 
those obtained by the Rayleigh-Riz method, but the former 
were somewhat smaller than the latter, which might be due 
to the effects of Poisson’s ratio. Therefore, this result sug¬ 
gests that all hygroexpansion forces are not stored in the 
strips but contribute to the deflection after buckling begins. 
In addition, this result indicates that the postbuckling be¬ 
havior of the strips differs greatly from that in the panel. 


Conclusions 

This study was carried out to investigate postbuckling be¬ 
havior of wood-based sandwich panels under various envi- 
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ronmental conditions. For the model development of 
postbuckling behavior, a few assumptions were made, in¬ 
cluding the free hygoexpansion that could be converted into 
elastic strain, the ignoring of imperfections, the time- 
dependent creep, and relaxation phenomena. Taking the 
edge movement of hardboard fixed on the core into consid¬ 
eration resulted in some discrepancy between theoretical 
and experimental results of the center deflection, and the 
theoretical value showed a slight overestimate because of 
neglecting the plastic deformation. This study showed that 
there was a large difference in the postbuckling behavior 
between panels and strips. The postbuckling behavior of 
panels depended on the MOE, while the behavior of strips 
was independent of MOE. Furthermore, the in-plane load 
of the panel increased beyond buckling, but the load for 
strips was constant with the initial buckling load. Further 
studies are necessary to improve the accuracy of the devel¬ 
oped model on buckling and postbuckling, including the 
viscoelasticity dealing with stress changes with time and the 
development of a more sophisticated model that accommo¬ 
dates the plastic deformation of panels and the nonuniform 
edge movements. 


Appendix: abbreviations and symbols 

n total potential energy of panel 

Ay, Dij extensional and bending stiffness, respectively 
e x , s y midplane strains in x and y directions, 

respectively 

y midplane shear strain in xy plane 

N x , N y hygroscopic forces in x and y directions, 

respectively 

h panel thickness 

a, b panel length in x and y directions, respectively 

x, y, z coordinate system 

u, v, w displacement in x, y, and z directions, 

respectively 

AMC moisture content change 

a„ ci y coefficient of linear expansion in x and y direc¬ 
tions, respectively 

LE linear expansion ( aAMC ) 

s arc length 

R radius of curvature 

N x , N y resultant in-plane forces per unit length in x and 

y directions, respectively 


N n nondimensional load, /Vf b 2 !E 2 ti 

N a nondimensional critical load, N^ a b 2 IE 2 h 3 

ft ratio of the core to the total width, aja + a c 

E c effective core MOE, E„ + E ± (i.e., the summa¬ 

tion of MOE parallel to the grain and perpen¬ 
dicular to the grain) 
h c core thickness 
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